NUWC-NPT  Technical  Report  10,275 
21  January  1993 


AD-A261  279 


Exact  Detection  Performance  of 
Normalizer  with  Multiple-Pulse 
Frequency-Shift-Keyed  Signals  in  a 
Partially-Correlated  Fading  Medium 
With  Generalized  Noncentral 
Chi-Squared  Statistics 


Albert  H.  Nuttall 
Surface  ASW  Directorate 


Naval  Undersea  Warfare  Center  Detachment 

New  London,  Connecticut 


Approved  for  public  release;  distribution  Is  unlimited. 


93  2  24  092 


PREFACE 


This  research  was  conducted  under  NUWC  Project  Number  A70272, 
Subproject  Number  RR00N00,  Selected  Statistical  Problems  in 
Acoustic  Signal  Processing,  Principal  Investigator  Dr.  Albert  H. 
Nuttall  (Code  302).  This  technical  report  was  prepared  with 
funds  provided  by  the  NUWC  In-House  Independent  Research  Program, 
sponsored  by  the  Office  of  Naval  Research. 


The  technical  reviewer  for  this  report  was  Keith  Peters  (Code 
3314  ) . 


REVIEWED  AND  APPROVED:  21  JANUARY  1993 


Mc¬ 
Donald  W.  Counsellor 
Director,  Surface  Antisubmarine  Warfare 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  HO  0704  0188 

I'  jOlif  rcourim'i  Dur'Mn  *»»r  thi%  •  Mmsfion  • )t  »nhjim.»fi>!P  <s  »<|p  1  G*1  '»*sp»>r\*  »<■  'uOinq  th*  l * rn*»  f . >r  M-vi#Ain  |  ti.'»n  -.r.if  hi»>  ]  i.»i  j  \..ur-  i»\ 

j.Hhfftft'l  in#t  fh*  'l-tt-i  n+fii+fj  tntj  ,  imf  n<  /inmntf  fhi*  ■!>#•*  **•.-#»  »>  »fi<  «  ,*ru)  •  r^< j  i  h.*  |ji;i  *  |rn  'Mim.it*  u  »n,  i»**i  t  t 

»f  mt-irmrtfi- *n.  .m  nj#linfj  '<)»  m-j  i*..-.  r*j  JV.iihmijt<jn  •  *»,.iOiju.»f  /h  rt  i  •«  *i  •  ir.itP  *•  ■>  m.inruii.m  ■  'i1*-  timns  imj  'i  .*■> 

j4vi%  Mi>jhM4y.  'V04  AMiivjIfin.  /A  iJH)J  A\l)}  ind  •<  <  tMi*  j*f«  *  it  M4n.i>j*m»«<t  infj  Hijdj**'  1*  »p^»  *tmk  Kr««jij<  Nun  »■"!,»-* !  iO/fM  0  tHH)  /V  ishtnqU ■»*  it 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  OATE  3.  REPORT  TYPE  AND  DATES  COWERED 

21  January  1993  Progress 

a.  title  and  subtitle  Exact  Detection  Performance  of 

Normalizer  with  Multiple-Pulse  FSK  Signals 
in  a  Partially-Correlated  Fading  Medium  with 
Generalized  Noncentral  Chi-Squared  Statistics 

S.  FUNDING  NUMBERS 

PE  61152N 

6.  AUTHOR(S) 

Albert  H.  Nuttall 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Naval  Undersea  Warfare  Center  Detachment 

New  London,  Connecticut  06320 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NUWC-NPT 

TR  10275 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

Chief  of  Naval  Research 

Office  of  Naval  Research 

Arlington,  VA  22217-5000 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release; 
distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

The  detection  and  false  alarm  probabilities,  for  a  normalizer 
operating  in  a  partially-correlated  fading  medium  with  an 
additive  unknown  noise  level,  are  derived  exactly  in  the  form  of 
series  expansions  which  can  be  quickly  and  accurately  computed. 

This  new  work  is  an  extension  of  [1],  where  the  characteristic 
function  of  the  decision  variable,  for  a  system  operating  in  a 
known  background  noise  level,  was  derived  exactly  and  evaluated 
for  a  variety  of  cases.  These  new  results  also  replace  the 
approximate  analysis  of  the  same  normalizer  system,  that  was 
conducted  in  [2],  with  an  exact  analysis  here. 

The  specific  situation  considered  here  is  as  follows.  The 
transmitted  signal  consists  of  K  pulses  separated  in  time, 
frequency  space  so  as  to  be  nonoverlapping.  In  passing  through 
the  medium  to  the  receiver,  each  signal  pulse  is  subjected  to 
fading.  In  particular,  pulse  pairs  which  are  closely  spaced  in 

14.  SUBJECT  TERMS 

detection  performance  normalizer 

multiple  pulses  frequency-shift-keyed 

partial-correlation  fading  medium 

15  NUMBER  OF  PAGES 

56 

16  PRICE  CODE 

17  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20  LIMITATION  Ci  ABSTRACT 

SAR 

**0-0 * -2 80-5500  *orm  *9l*  *»w 


■'s  '  r*>U  P>  -•  ‘4  V  -t  1  i  •* 


;  *<*' 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION 

OF  THIS  PAGE _ _ _ 

13.  ABSTRACT  (continued) 

time,  frequency  space  can  fade  in  a  highly  dependent  fashion, 
while  those  pairs  more  widely  separated  can  have  relatively 
independent  fading  behavior;  that  is,  the  transmitted 
frequency-shift-keyed  signal  pulses  undergo  partially- 
correlated  fading  of  a  very  general  character  that  can  contain 
both  deterministic  components  as  well  as  random  components. 
Thus,  for  example,  the  amplitude- fading  statistics  are  not 
limited  to  be  Rayleigh. 

Additive  zero-mean  Gaussian  noise  of  unknown  level,  which 
is  stationary  over  the  total  signal  transmission  time  and 
which  has  a  flat  spectrum  over  the  total  signal  bandwidth,  is 
present  at  the  input  to  the  receiver,  in  addition  to  the 
fluctuating  signal  pulses  (when  present). 

Receiver  processing  consists  of  matched  filtering  of  each 
of  the  K  time-delayed,  frequency-shifted  pulse  locations, 
followed  by  squared-envelope  detection,  sampling,  summation, 
and  comparison  of  this  decision  variable  with  a  variable 
threshold,  for  statement  on  signal  presence  or  absence.  The 
variable  threshold  is  obtained  by  estimating  the  noise  power 
in  L  time,  frequency  bins  in  which  signal  is  known  to  be 
absent.  The  additional  uncertainty,  caused  by  this  lack  of 
knowledge  of  the  additive  noise  level,  requires  additional 
input  signal-to-noise  ratio  to  maintain  the  same  performance 
as  for  the  f ixed-threshold  (known  noise  level)  system.  The 
exact  quantitative  evaluation  of  this  degradation  is  the 
subject  of  this  study. 
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EXACT  DETECTION  PERFORMANCE  OF  NORMALIZER  WITH  MULTIPLE-PULSE 
FREQUENCY-SHIFT-KEYED  SIGNALS  IN  A  PARTIALLY-CORRELATED  FADING 

MEDIUM  WITH  GENERALIZED  NONCENTRAL  CHI-SQUARED  STAT!  TICS 

INTRODUCTION 

In  a  recent  study  [1],  the  detection  performance  capability 
of  a  multiple-pulse,  frequency-shift-keyed  system,  with  received 
signals  subject  to  partially-correlated  fading,  was  analyzed 
exactly  in  terms  of  the  characteristic  function  of  the  system 
output  decision  variable.  It  was  assumed  there  that  the  additive 
background  noise  level  was  known,  so  that  a  fixed  threshold  could 
be  set  for  an  arbitrarily  specified  false  alarm  probability. 

Then,  the  detection  probability  was  evaluated  exactly,  as  a 
function  of  the  threshold  level,  the  received  signal-to-noise 
ratio,  the  number  K  of  signal  pulses,  and  the  fading  statistics. 

Here,  we  extend  these  earlier  results  to  cover  the  case 
where,  instead,  the  background  noise  level  is  unknown  and  must  be 
estimated  on  the  basis  of  a  finite  number  L  of  noise-only 
samples.  Furthermore,  the  current  analysis  will  also  be  exact, 
yielding  new  series  expansions  for  the  detection  and  false  alarm 
probabilities.  These  new  results  replace  the  approximate 
analysis  of  the  same  normalizer  system  that  was  conducted  in  [2]. 
This  normalizer  possesses  the  important  capability  of  constant 
false  alarm  rate;  that  is,  a  prescribed  false  alarm  probability 
can  be  achieved,  exactly,  in  the  absence  of  knowledge  of  the 
actual  noise  background  level. 
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The  reader  is  referred  to  [1]  for  additional  background, 
motivation,  interpretations,  and  related  references.  For  the 
sake  of  brevity,  we  will  employ  the  same  notation  as  in  [1]  and 
presume  that  the  reader  has  complete  familiarity  with  the  earlier 
material  and  developments. 
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PROBLEM  DEFINITION 


BASIC  SIGNAL  PROCESSING 

The  problem  of  interest  will  be  couched  in  a  time,  frequency 
setting,  in  which  the  individual  transmitted  signal  pulses  occupy 
disjoint  nonoverlapping  cells;  however,  the  extension  to  other 
variables,  such  as  space  or  angle,  should  be  obvious. 

A  set  of  K  signal  pulses  is  transmitted  at  known  time, 
frequency  locations,  as  indicated  in  figure  1.  These  pulses  may 
be  abutting  in  time  and/or  frequency,  or  they  may  be  widely 
separated  in  time,  frequency  space.  The  individual  signal  pulses 
need  not  have  unity  time-bandwidth  product,  although  many 
practical  applications  would  utilize  that  format. 

frequency 


Figure  1.  Time,  Frequency  Occupancy  Diagram 


3 


TR  10275 


At  the  receiver,  the  noisy  received  waveform  is  processed 
through  K  matched  filters,  which  utilize  knowledge  of  the 
particular  time,  frequency  occupancy  pattern  of  the  signal 
pulses.  Their  squared  envelope  outputs  are  then  appropriately 
sampled  and  summed  for  comparison  with  a  scaled  estimate  of  the 
background  level  obtained  from  the  reference  collection  of  L 
roise-only  bins.  The  scale  factor  applied  to  the  noise-level 
estimate  yields  the  threshold  value,  which  is  selected  to  realize 
the  desired  false  alarm  probability.  If  the  threshold  is 
exceeded  by  the  sum  of  K  squared  envelopes,  a  signal  is  declared 
present  in  that  pattern  of  K  pulses. 

Depending  on  the  time,  frequency  separation  of  the  signal 
pulses  in  figure  1,  the  received  signal  strengths  of  the 
individual  pulses,  when  signal  is  present,  may  fade  considerably 
and  with  statistical  dependence  between  them.  The  exact  amount 
and  effect  of  the  fading  depends  on  the  distribution  of  the 
fading  and  on  the  Kxk  covariance  matrix  between  the  fading 
amplitudes  applied  to  each  signal  pulse. 

It  is  presumed,  as  in  [1],  that  during  the  time  extent  of  a 
single  pulse,  the  fading  of  that  particular  pulse  is  essentially 
constant,  resulting  in  a  constant  amplitude  scaling  and  phase 
shift  applied  to  each  received  pulse.  The  phase  shifts  of 
distinct  pulses  can  have  arbitrary  statistics  and  dependencies; 
the  particular  phase  shifts  are  irrelevant  in  the  method  of 
signal  processing  adopted  here  and  in  [1]. 
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The  choice  of  L,  the  number  of  noise-only  pulses,  must  be 
based  on  a  compromise.  Larger  L  results  in  a  more  stable 
estimate  of  the  background  noise  level,  in  which  case  the 
performance  of  the  normal izer  processor  for  figure  1  approaches 
that  for  the  fixed-threshold  processor  in  [1].  However,  L  cannot 
be  chosen  so  large  as  to  encounter  nonstationary  and/or  nonwhite 
noise  fields  over  the  time, frequency  extent  depicted  in  figure  1. 

For  signal  pulses  consisting  of  bursts  of  pure  tones,  the 
receiver  implementation  of  figure  1  can  be  accomplished  by  means 
of  fast  Fourier  transforms.  This  case  would  utilize  signal 
bursts  and  noise  slots  that  have  approximately  unity  time- 
bandwidth  products,  and  would  result  in  minimum  consumption  of 
time,  frequency  space  in  figure  1.  However,  the  results 
presented  below  are  not  limited  to  this  special  case. 
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REVIEW  OF  PARAMETERS 

Very  heavy  reliance  will  be  made  upon  the  methods  and  results 
presented  for  the  f ixed-threshold  system  analyzed  in  [1],  and 
which  the  reader  must  be  familiar  with.  However,  it  is  very 
worthwhile  at  this  point  to  review  and  summarize  the  most 
important  and  critical  parameters  that  governed  performance  in 
[1],  and  which  are  relevant  to  the  normalizer  analysis  here,  as 
well.  These  parameters  are  listed  below;  their  first  appearance 
and  detailed  explanation  can  be  found  in  [1;  pages  iv  -  viii]  in 
the  list  of  symbols  there. 


lm 


o 

(m) 


(m) 


x(m) 

-k 


D 


km 


,(m) 

k 

r(m) 

kj 

.(m) 


m 


number  of  transmitted  signal  pulses;  see  figure  1  above 

number  of  components  in  fading  model;  see  [1;  (45)  -  (47)] 

average  random  received  signal  energy  in  the  m-th 
component  of  a  single  pulse,  1  <  m  £  M;  see  [1;  (52)] 

one-sided  spectral  level  of  received  white  noise  (unknown) 

Kxk  normalized  covariance  matrix  of  signal  amplitude 
fading;  see  [1;  (127)] 

kxk  eigenvalue  matrix  of 

k-th  eigenvalue  of  matrix  A^m^,  1  <;  k  <  K 

deterministic  received  signal  energy  in  the 
m-th  component  of  the  k-th  pulse;  see  [1;  (50)] 

k-th  eigenvector  of 

j-th  component  of  vector  V^m^,  1  £  j  <  K 
auxiliary  constants;  see  [1;  (137)] 

fractional  strengths  of  random  components;  see  [1;  (191)] 
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It  is  also  worthwhile  at  this  point  to  expound  on  the 
covariances  of  the  fading  model,  since  these  will  be  used  at  a 
later  stage  of  the  development.  The  covariance  function  of  the 
alternating  component  gm(t,f)  in  [1;  (45)  &  (47)]  is 

Rmn<'r'v>  =  Vt,f)  9n(t-T,f-v>  ,  (1) 

while  the  covariance  of  power-scaling  variate  qR  =  is, 

from  [1;  (43),  (44),  (54),  (55)] 

Rkj  =  Cov(qR,qj)  =  Cov [q( tR, f R) ,  q( t ^ , f ^ ) )  .  (2) 

The  covariance  coefficient  between  qR  and  q^  is 


from  [1;  (56),  (66)  or  (69)].  We  also  define  the  normalized 
covariance  function 


COV ( T , V ) 


Rn(T,^) 

Rn(0,0) 


(4) 


as  in  [1;  (194)],  from  which  there  follows 

Pkj  -  (5) 

as  in  [1;  (56),  (66)  or  (69)].  Finally,  we  define  a  special  case 
of  (4),  for  the  case  of  equally  time-spaced  signal  pulses  without 
any  frequency  shifts,  according  to 
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Rll( tk+l"tv'°) 


covj  -  c^k-n-V0*  =  — r'1  ~ ( 0T0) 


for  all  k 


(6) 


since  this 
in  ( 6 )  was 


quantity  will  be  very  useful  later;  this  last  quantity 
denoted  by  Cov  in  [1]  and  all  the  figures  therein. 
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EXCEEDANCE  DISTRIBUTION  FUNCTION  OF  NORMALIZER  OUTPUT 


FIXED-THRESHOLD  PROCESSOR  OUTPUT 


Before  we  begin  the  analysis  of  the  normalizer,  we  need  to 
review  and  manipulate  the  results  for  the  fixed-threshold  case. 
The  general  result  for  the  characteristic  function  of  signal 
processor  output  y,  in  the  case  of  M  uncorrelated  components  in 
fading  model  [1;  (45)],  is  given  by  [1;  (136)]  as 


£y(M 


-  (1  -  i2t)K(i,M-11 


exp 


M  K 


m=l  k=l 


where 


2E. 


0  =  i  +  _ IB  x(m) 

ymk  1  N  -k 
o 


emk  =  * 


M 

K  / 

TT 

m _  1 

n 1  ■ 

i2t2mk) 

k  m-  j. 

emk  1 

i  - 

i2£^kJ 

9 

2 

(m) 

< 

for  1 

1  m  <  M 

(7) 


(8) 


This  situation  pertains  to  ignoring  the  L  noise-only  outputs  in 
the  lower  portion  of  figure  1. 

We  will  limit  consideration  here  to  the  special  case  6 
delineated  in  [1;  page  23],  Then,  eigenvalues  (X^J  and 
eigenvectors  {V^j  are  not  functions  of  component  number  m,  and 


ymk  "  1  + 


2E1 

N  ^m  — k 
o 


from  (8)  and  [1;  (191)].  Also,  from  [1;  (157)], 


(9) 


9 


TR  10275 


K 


(10) 


We  now  want  to  perforin  a  transformation  of  variable  l  in 
characteristic  function  (7).  In  particular,  consider  the  use  of 
the  slightly  more  general  transformation  (relative  to  that  used 
in  [1;  (D-3)]),  where  scale  factor  p  is  under  our  choice: 


z 


(P  >  0)  . 


(11) 


Substitution  into  (7),  and  expansion  in  a  power  series  in  z 
(according  to  appendix  A),  yields  the  characteristic  function  of 
the  fixed- threshold  processor  output  y  in  the  (modified  but 
exact)  form 


f  Al)  =  F  zK  fn  zP  “  F  III  fn  (1  ~  i^)"K'P  '  (12> 


p=0 


p=0 


where  the  constants  are  found  according  to 


=  {§.} 

[2) 


K 


M  K 

TIT 7  Q 

m=l  k=l 


mk 


a_  -  -  o 


!  M  K  emk 


2  L _ ■  L - -  o  ' 

m=l  k=l  ymk 


(13) 


M  K 


P  a. 


airr  [i.  -i—i 

2  m=l  k=l  '  2^mk 


P-1 


_  _A 


2Q 


mk 


1  ~  P 


'mk 

^mk 


-  k  ("  -  i)  t1  -  !)P  £or  p  2  i  - 


(14) 


fo  *  exP<ao> 


f  =  -  )  n  a  f  for  p  1  1 

P  P  fjrf  n  p-n 


(15) 
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For  p  =  2,  (14)  reduces  to  [1;  (185)  or  ( D— 10)]-  However,  we 
want  to  choose  g  here  so  that  the  magnitudes  of  the  coefficients 
(F  f  }  in  (12)  decay  as  rapidly  as  possible  with  p.  On  the  other 

Jr 

hand,  some  of  the  { f ^ }  can  then  be  negative. 

Now,  suppose  we  have  found  the  characteristic  function  of  any 
random  variable  y  in  the  series  form  (12),  where  K  is  integer. 
Then,  the  origin  value  of  the  characteristic  function  is 
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E  (u)  =  Pr( y  >  u)  .  (20) 

Here,  in  (18),  we  used  the  following  definition  of  the  exceedance 
distribution  function  E,  which  can  be  expressed  in  terms  of 
incomplete  gamma  functions: 


E(v,u)  = 


v-1 

x  exp 
F(v) 


r(  v ) 


(21) 


More  details  on  recursions  for  these  functions,  when  v  is  equal 
to  an  integer  or  a  half-integer,  are  given  in  appendix  B. 
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NORMALIZER  OUTPUT 

Now,  let  fixed  threshold  u  in  (20)  be  replaced  by  a  variable 
threshold  according  to 

u  =  v  yq  ,  (22) 

where  v  is  a  positive  scaling  constant,  and  yQ  is  an  independent 
positive  random  variable  given  by  the  sum  of  the  L  independent 
noise-only  bin  outputs  in  figure  1.  This  latter  quantity  is 
proportional  to  an  estimate  of  the  noise  background  power  level 
in  which  the  normalizer  system  is  operating.  The  characteristic 
function  of  yQ  can  be  found  from  (1;  (40)],  by  replacing  K  by  L 
and  by  setting  signal  strength  E  to  zero: 

f  (1)  =  (1  -  i2l)~h  .  (23) 

ro 

The  corresponding  probability  density  function  of  yQ  is 

pv  (x)  =  - -  expJ_-x/2 )  for  x  >  o  .  (24) 

Yo  (L-l) i 

We  are  now  interested  in  the  exceedance  distribution  function 
En(v)  of  the  normalizer  output,  namely  ratio  y/yq: 
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F  ZZ  fn  f  dx  xL~t1  exP<- 
p=0  P  l  2L  ( L— 1 ) 


00  K-l+p  ,  ,  *k  r  L-l+k  ,  , , 

=  F  ZZ  fD  C  |r  I  fix  -i -  exp  -X  i  +  |  = 

p=0  P  k=0  k*  {fiJ  {  2L  ( L- 1 )  !  1  12  P}) 


L-l+k 


2  (L-l) 


£ _ ]L  f  Ky^P  [L-l+k|  r 2v 1K 

+  2vJ  £p  far  l  *  )  Im  2vJ 


L  00 


F  (irb*)  £p  T(K-in..L-i,r|3 


(25) 


where  we  define  the  auxiliary  functions 


k  in 


HI  /  .  1,  -I  i,  m  „J\.  HI 

T(m,n,x)  h  YIZ  k  «  =  IZ  (n+l)k  fy  =  IZ  U(k,n 

k=0  v  K  '  k=0  K  K*  k=0 


/X)  .  (26) 


As  a  partial  check  on  (25),  we  have  origin  value 


E_(0)  =  F  2 _  fn  T ( K- 1+p, L-l , 0 )  =  F  £ _ ’  f  =  1  , 

p=0  P  p=0  P 


(27) 


using  ( 16 ) . 


For  rapid  evaluation  of  functions  T  and  U  in  (26),  we  have 


recursions 


T(m,n,x)  =  T(m-l,n,x)  +  U(m,n,x)  for  mil. 


(28) 


U(m,n,x)  =  U(m-l,n,x)  x  for  mil, 


(29) 


along  with  starting  values 


T(0,n,x)  =  1  ,  U(0,n,x)  =  1 


(30) 
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A  general  result  that  will  find  frequent  application  is 


( l~x)n+1  T(m,n,x)  =  1  -  xm+^  T(n,m,l-x)  . 


(31) 


As  a  special  case,  we  also  have  [3;  15.1.8] 

oo  k 


T(°°,n,x)  =  2 _  (n+l)R  FT  =  iFo^n+1;x)  =  ( 1  ~  x> 

k=  0 


-n-1 


(32) 


By  means  of  this  last  relation,  we  have 


f-Ll  2v 

1  ?v  1 

-L 

[3  +  2vl 

L 

i  ,L  1,(3  +  2vJ 

r  p  +  2vJ 

P 

/ 

(33) 


which  enables  the  normalizer  exceedance  distribution  function  in 
(25)  to  be  expressed  in  the  alternative  form 


En<v>  - F  (<£dL  n:  £p  -  (er2)L  +  (‘x2)1]' 


-  1  -  F  fp  [l  -  (gf^)  t(k-1+P,L-1/?|^)]  ,  (34) 

where  we  used  (16).  This  latter  form  is  recommended  over  (25) 
because  both  f^  and  the  bracketed  term,  [  ],  in  (34)  go  to  zero 
as  p  ->  °°.  The  approach  of  f  to  zero  follows  from  (16).  It  is 
also  very  important  to  notice  that  the  bracketed  term  in  (34)  is 
not  a  function  of  the  signal-to-noise  ratio. 

By  means  of  (31),  (34)  can  also  be  expressed  as 


00  I  9  \K+P  (  ft  x 

En(v)  -  1  -  F  IZ  fp  (j^)  T  (l_  1  ,E“  1+P' p+2v)  •  <35> 
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Although  both  f  and  the  succeeding  power  in  (35)  do  tend  to  zero 
as  p  ■»  ®,  function  T  does  not;  of  course,  the  product  of  the 
power  term  and  T  do  tend  to  zero.  Nevertheless,  (34)  will  be 
more  useful,  since  a  simple  recursion  on  p  exists.  (It  has  been 
numerically  verified  that  (25),  (34),  and  (35)  are  all 
independent  of  the  particular  value  of  the  positive  scaling 
constant  p  in  transformation  (11);  we  are  now  free  to  choose  the 
value  of  p  that  results  in  the  most  rapid  convergence  of  the 
various  series.)  Results  (25),  (34),  and  (35)  are  alternative 
expressions  for  the  detection  probability  PD(v)  of  the  normalizer 
output  ratio  Y^Yqi  with  threshold  value  v;  that  is, 

PD( v)  =  Pr(y/Y0  >  v)  for  fixed  v. 

The  numerical  evaluation  of  form  (34)  requires  a  difference 
in  the  bracketed  term  and  thereby  suffers  a  loss  of  accuracy.  A 
more  useful  numerical  alternative,  using  a  sum  of  positive  terms, 
is  furnished  by 


(36) 


This  result  follows  upon  use  of  (26)  and  (33).  The  inner  sum  on 

k  in  (36)  obviously  goes  to  zero  as  p  increases,  as  does 

coefficient  f  .  Recursions  for  the  inner  sum  on  k  make  this  a 
P 

rapid  accurate  alternative  to  (34);  result  (36)  is  programmed  and 
used  in  this  investigation. 
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CHECK  ON  EXCEEDANCE  DISTRIBUTION  (25) 


Result  (25)  gave  the  normalizer  exceedance  distribution 
function 

En(v)  =  Pr( y  >  v  yo)  *  Pr(^“  >  v)  *  (37) 

If  we  now  let  the  normalizer  fixed  threshold  v  be  given  by 


v 


(38) 


then  we  find 


En(v) 


=  Pr(Y  >  u  I?)  * 


(39) 


But  as  L  ->  ®,  random  variable  yq/(2L)  tends  to  the  value  1,  since 
the  mean  of  yq  is  2L  and  its  standard  deviation  is  2/L;  see  (24), 

Ir 

for  which  the  k-th  cumulant  is  L  (k-1)!  2  .  Thus,  we  should  have 


lim  E 
L->°° 


n 


Pr ( Y  >  u)  =  Ey(u)  , 


(40) 


where  the  latter  exceedance  distribution  was  given  by  (18). 

In  order  to  confirm  this  result,  substitute  v  =  u/(2L)  into 
(25)  to  obtain 


En(5?'  * 


;)  * F  (r^)"  §  fp  t(k-1+P'l-1'p^7l)  • 


(41) 


But 


(s“A/l) 


(r 


+  u/((3L) 


t 


as  L  -»  ® 


(42) 
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while 


(r^)k  = 


m  ( L) ,  ,  /0 .  k  m  ,  /..k 

r - -  J _ Me  (u/g  )  _  v - -  (  u/g  ) 

L .  V  i  b  ~  L .  k  i 

k=0  K'  (L  +  u/0)K  k=0  K* 


as  L  -»  °°  . 


(43) 


Therefore , 


lim  E  [^7] 
L-*»  nULj 


=  F  exp 


in  agreement  with  (18).  Here,  we  used  (B-17). 


00 

K-l+P 

EZ 

fp  ^  ' 

p=0 

p  k=0 

Here 

,  we  used 

k^  *  F  fD 

K*  p=0  p 


E(K+P,a) 

(44) 


18 


TR  10275 


FALSE  ALARM  PROBABILITY 


If  signal  is  absent  from  the  processor  input,  the  earlier 
results  for  the  normalizer  output  exceedance  distribution 
function  should  reduce  to  the  false  alarm  probability.  The 
moments  of  normalizer  output  y/Y0  ate  presented  in  appendix  C, 
when  signal  is  absent.  Also,  in  this  case,  we  have  Qmk  =  0, 


'mk 


=  0,  giving,  from  (13)  -  (14), 


(  p 

F  =  £  a  = 

*  \2)  '  o 

0  , 

p 

a  =  K 
P 

(‘  - 1) 

for  p  1  1 

There  follows  (see  (A-5) 

-  ( A-6 ) ) 

00 

C  fD  zP  =  exp 

p=0  p 

00 

ho 

a  zp 
P 

=  exp 

K  T~~*  — 

p=t  p 

(>  -  f)P  ‘p 

=  exp |-K  In | 1  - 

i 

Ml- 

b-(‘- 

f)z]  K  - 

(45) 


from  which  we  find  coefficients 


(K) 
f  =  — 
P  P 


'd 

7^  1  -  §  for  p  1  0  . 


(47 


Then  (25)  reduces  to  the  normalizer  false  alarm  probability 


pf<v>  -  (f)K  (r^)L  -  if  T(K-1+p-L-1'rTTv) 


p=0 


(48) 
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Furthermore,  since  this  result  must  be  independent  of  p,  we  set 
P  =  2  to  get 


PF(v)  = 


(1  +  v) 


T  (K-1',L-1/T—+— v)  - 


{ti  <L>k 


(1  +  v)“  k=0 


k! 


,  \k 

rH  - 


(49) 


=  1 


-  h-N 


K  L-l  (K). 


ji 


j=0 


(1  +  v) 


(50) 


where  we  used  (31).  This  result,  (50),  agrees  with  [2;  (22)]. 

The  form  in  (49)  utilizes  a  sum  of  all-positive  terms  and  thereby 
retains  significance  for  very  small  Pf. 
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GRAPHICAL  RESULTS 


The  scenario  to  be  considered  here  is  the  same  as  that  used 
in  [1]/  namely,  a  sequence  of  K  signal  pulses,  equally-spaced  in 
time  but  with  no  frequency  shifts.  Also,  the  covariance  between 
the  signal  amplitude  fading  coefficients  will  be  taken  to  be 
exponential  in  time  separation.  In  this  manner,  we  can  make  a 
direct  comparison  of  the  results  for  the  fixed-threshold 
processor  in  [1]  with  the  normalizer  here,  and  numerically 
determine  the  degradation  associated  with  utilizing  too  small  a 
value  for  L,  the  number  of  noise-only  bins.  As  in  [1;  page  68], 
the  KM  deterministic  signal-to-noise  ratios  are  all 

taken  to  be  zero,  and  the  M  strengths  { r ^ m ^ }  are  all  1,  for  the 
numerical  examples  investigated  here. 

Due  to  the  multitude  of  parameters  in  this  investigation,  it 
is  expedient  to  hold  all  but  one  of  them  fixed,  and  to  let  one 
parameter  vary,  thereby  determining  the  effect  of  that  parameter 
by  itself.  The  first  series  of  plots  in  figures  2,  3,  4,  5,  6 
correspond  to  varying  K,  the  number  of  signal  pulses  in  figure  1, 
over  the  range  K=2,  3,  4,  5,  6.  The  curve  for  K  =  1  has 
already  been  given  in  [2;  figure  32];  that  analysis  was,  in  fact, 
exact  for  K  =  1,  M  =  2. 


The  abscissa  in  all  the  plots  is  L,  the  number  of  noise-only 
bins  that  are  used  to  estimate  the  background  noise  level,  and 
which  varies  over  the  range  from  1  to  100  in  all  cases.  The 
asymptotic  values  of  required  signal-to-noise  ratio  at  L  =  ®  are 
available  in  [1].  The  ordinate  in  the  following  plots  is  the 
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input  signal-to-noise  ratio  (SNR)  E1/Nq  in  dB,  required  in  order 
to  realize  the  specified  performance  levels  in  terms  of  detection 
and  false  alarm  probabilities,  where  Ej  is  the  average  random 
received  signal  energy  in  one  pulse,  and  Nq  is  the  additive  noise 
power  density  level.  The  curves  are  parameterized  in  terms  of 

O  A 

the  false  alarm  probability,  which  is  set  to  values  10”  ,  10~  , 

10  ,  and  10“  .  The  particular  parameter  values  that  apply  to 

each  plot  appear  in  the  title  for  each  figure.  Parameter  Cov^ 
was  defined  and  explained  in  (1)  -  (6);  basically,  it  is  the 
covariance  between  adjacent  signal  pulses. 

For  most  rapid  convergence  of  the  series  for  the  exceedance 
distribution  function  of  the  normalizer  output,  the  value  of 
scaling  p  in  transformation  (11)  should  be  changed  for  each  point 
computed  on  these  plots.  Some  samples  for  the  best  p  values  are 
tabulated  in  appendix  D,  along  with  the  particular  program  used 
for  computing  the  results  in  all  the  figures  here. 

The  first  result  in  figure  2,  for  K  =  2,  replaces  the 
approximate  result  presented  in  (2;  figure  33].  It  will  be  seen 
that  approximately  1  dB  less  is  required,  according  to  the  exact 
results  of  figure  2  here;  thus,  the  approximate  analysis  in  [2] 
was  somewhat  pessimistic  in  terms  of  the  input  signal-to-noise 
ratio  required  to  realize  a  specified  performance  level.  This 
same  conclusion  was  also  reached  upon  a  comparison  of  the  exact 
analysis  of  the  fixed-threshold  processor  in  [1]  with  the 
corresponding  approximate  analysis  in  [4]. 
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As  the  number  of  pulses,  K,  increases  from  2  to  6  in  figures 
2  through  6,  the  required  value  of  per-pulse  SNR,  E^/Nq, 
decreases  monotonically .  However,  the  total  SNR,  summed  over  all 
K  pulses,  increases  with  K;  this  same  effect  was  observed  and 
investigated  quantitatively  in  [4],  although  by  an  approximate 
analysis. 

The  next  series  of  plots,  in  figures  7  -  11,  correspond  to 
varying  covariance  coefficient  Cov^^  in  (6)  over  the  range  0  to  1 

i* 

in  steps  of  (n/4)  ,  for  n=0,  1,  2,  3,  4.  The  remaining 
parameters  are  held  fixed  at  values  PD=  .5,  K=4,  M=2.  As 
expected,  the  required  input  signal-to-noise  ratio  increases 
monotonically  with  Covj.  The  increases,  as  n  makes  the  four 
changes  from  0  to  4,  are  in  the  ranges  .15  to  .21  dB,  .21  to  .28 
dB,  .31  to  .35  dB,  and  .36  to  .45  dB,  respectively.  The  total 
increases,  as  Covj  varies  from  0  (figure  7 )  to  1  (figure  11)  are 
in  the  range  1.1  to  1.2  dB.  Thus,  at  least  for  this  example,  the 
cost  of  the  fading  covariance  coefficient  approaching  unity  is 
moderate . 

The  final  series  of  plots,  in  figures  12  -  14,  depict  the 
variation  of  required  input  signal-to-noise  ratio  with  M,  the 
number  of  components  in  the  fading  model  [1;  (45)  -  (47)]. 

(The  plot  in  figure  12  is  identical  to  figure  4  above.)  As  M 
changes  from  2  to  3,  with  other  parameters  fixed,  the  required 
input  SNR  decreases  by  approximately  1  dB,  while  the  change  of  M 
from  3  to  4  allows  for  an  additional  decrease  of  . 5  dB  in 
required  input  SNR.  Generally,  larger  M  means  a  fading  channel 
model  subject  to  less  severe  deep  fades;  see  [1;  figures  2-8]. 
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Figure  3.  Required  SNR  for  PD  =  .9,  K=3,  M=2,  Cov^  =  /.5 


24 


TR  10275 


Figure  5.  Required  SNR  for  PD  =  .9,  K  =  5,  M  =  2,  Cov1  =  /.5 
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Figure  8.  Required  SNR  for  PQ  =  .5,  K=4,  M  =  2,  CoVj  *  /.25 


Figure  9.  Required  SNR  for  Pp  *  .5,  K  =  4,  M  =  2,  Cov^  =  / . 5 
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Figure  10.  Required  SNR  for  P 


5 ,  K 


2 ,  Cov 
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SUMMARY 

The  false  alarm  and  detection  probabilities,  for  a  normalizer 
operating  in  a  fading  medium,  have  been  derived  exactly  in  the 
form  of  series  expansions  which  are  capable  of  efficient 
recursive  numerical  evaluation.  These  results  allow  for  multiple 
signal  pulses  with  an  arbitrary  covariance  matrix  between  the 
amplitude  fadings  applied  to  each  signal  pulse.  Lack  of 
knowledge  of  the  additive  background  noise  level  requires  the  use 
of  a  group  of  noise-only  bins,  from  which  an  estimate  of  the 
noise  level  is  made  and  used  to  establish  a  threshold  for  the 
normalizer,  in  order  to  realize  a  specified  false  alarm 
probability.  The  degradation  associated  with  this  noisy 
threshold  has  been  evaluated  quantitatively;  a  program  for  its 
evaluation  in  a  specific  environment  is  furnished  in  appendix  D. 
Some  examples  of  its  use  have  been  presented  in  the  figures 
herein. 

The  expansion  of  the  characteristic  function  in  a  series,  as 
developed  here  in  appendix  A,  will  be  extended  in  [6]  to  include 
a  more  general  form  of  transformation  than  (11).  This  should 
allow  for  series  with  more  rapid  convergence  properties  than 
those  derived  here  in  (34)  -  (36),  for  example.  In  particular, 
extensions  of  the  Hermite  series  and  the  generalized  Laguerre 
series  used  in  [5]  will  be  developed. 
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APPENDIX  A.  SERIES  EXPANSION  OF  CHARACTERISTIC  FUNCTION  (7) 


The  transformation  of  variable  in  (11)  is  given  by 


1  4  r  _  1  -  z 

z  1  -  H0  '  ^  fiz  ' 


( A- 1  ) 


which  yields 


1  -  i2K  = 


zd  -  hi 

h&z 


1  -  i2lQ  = 


_  1  -  2(1  -  W( 


hfiz/Q 


.  ( A-2 ) 


When  these  expressions  are  substituted  into  characteristic 
function  (7)  and  simplified,  there  follows 


fY(M  =  F  z1 


1  -  z 


i  _  wn^lW-D 


M  K 


rrrrl1  -  -(i  -  nj] 

m=l  k=l 1 


x  exp  -  ( 1  -  z ) 


=1  k-1  ~  z{!  -  WOmk)  J  ' 


(A-3) 


where  F  is  given  by  (13). 

Except  for  the  factor  F  z  ,  we  now  take  the  natural  logarithm 
of  the  right-hand  side  of  (A-3)  and  expand  in  a  power  series  in 
z.  After  considerable  simplification  and  grouping  of  similar 
terms  together,  we  find  that  this  logarithm  is  given  by 


ap  Z 


(A-4) 


where  the  coefficients  { a^ 1  are  given  in  (13)  and  (14). 

Recalling  the  factor  F  z  in  (A-3),  we  obtain  the  characteristic 
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function  of  the  fixed-threshold  processor  output  y  in  the  form 

00 

f  (£)  =  F  zK  exp  YIZ  an  zP  •  (A-5) 

Y  lp=0  p 

We  now  employ  [5;  page  93]  to  obtain  the  desired  series  in 
powers  of  z  according  to 

£v(t)  <*-«> 

Y  p=0  p 

where  z  is  given  by  (A-l),  and 

f„  =  exp(a  )  ,  f  =  ^  >  .  n  a  f  for  p  >  1  .  (A-7) 

o  ^  o'  P  P  n  p-n 

Expression  (A-6)  is  the  desired  result.  Numerous  recursions 
can  be  developed  from  these  equations  that  enable  rapid  efficient 
numerical  evaluation  of  characteristic  function  f^.  More  will  be 
said  on  this  topic  in  the  program  for  evaluation  of  the  detection 
probability  for  the  normal izer  output;  see  appendix  D. 
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APPENDIX  B.  INCOMPLETE  GAMMA  RELATIONS 

Define  the  function 


C(v,x)  =  J  dy 


V—  1 

y  exp(-y) 
P(v) 


for  x  >  0  ,  v  >  0 


(B-l) 


(This  is  P(v,x)  in  [3;  6.5.1].)  Then,  from  [3;  6.5.2  and  6.5.3], 


C(v  x)  =  giy-jg  =  i  _  rtV.<-*I 

r( v)  1  r(v )  ' 


( B-2  ) 


in  terms  of  incomplete  gamma  functions.  The  function  C(v,x)  is  a 
cumulative  distribution  function;  that  is, 


C(v,«)  =  l 


( B-3 ) 


It  is  also  useful  to  notice  that 


C ( °°, x )  =  0  ,  ( B-4  ) 

since  the  unit  area  of  the  integrand  of  (B-l)  moves  progressively 
farther  out  in  the  positive  direction  as  v  increases.  This  can 
be  seen  from  the  fact  that  the  density  (integrand)  in  (B-l)  has 
mean  v  and  standard  deviation  /v. 

Also,  define  function 

\>- 1  .  _  . 

p(v,x)  =  - - r^)  ~X  for  x  >  0  ,  v  >  0  .  (B-5) 

Then  function  (B-l)  satisfies  recurrence 
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x 

C(v,x)  =  |  dy  p(v,y )  =  C(v-l,x)  -  p(v,x)  for  v  >  1  , 
0 

while 


p{v,x)  =  p(v-l,x)  for  v  >  1  . 

Convenient  starting  values  for  (B-6)  and  (B-7)  are 


( B-6 ) 


(B-7) 


C(Jj,x)  =  2  ♦  ( ( 2x ) **)  -  1  ,  C(  1  / x )  =  1  -  exp(-x)  , 

_L 

p (h,x)  =  (nx)  ^  exp(-x)  ,  p(l,x)  =  exp(-x)  . 

As  an  example,  with  v  =  n, 

p(n,x)  =  p(n-l,x)  for  n  =  2,3,4, — 


( B-8 ) 
( B-9 ) 


( B-10 ) 


C ( n, x )  =  C(n-l,x)  -  p(n,x)  for  n  =  2,3,4,...  (B-ll) 

along  with  the  starting  values  for  p(l,x)  and  C(l,x)  given  in 
(B-9)  and  (B-8),  respectively. 

As  a  second  example,  with  v  =  n  + 


p(n+*},x )  =  p(n-J*,x)  for  n  =  1,2,3, —  (B-12) 

Cfn+^x)  =  C(n-Jj,x)  -  p(n+1s,x)  for  n  =  1,2,3,...  (B-13) 

along  with  the  starting  values  for  p(h,x)  and  C(*s,x)  given  in 
(B-9)  and  (B-8),  respectively. 
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We  also  define  exceedance  distribution  function 


E(v,x) 


=  1  -  C(M,X)  =  J 


dy 


y^'1  exp(-y) 
T(v) 


=  r(v,x)/r(\>)  for  x  >  0  ,  v  >  0 


Then 


E(v,x)  =  E(v-l,x)  +  p(v,x)  for  v  >  1  , 


with  recursion  (B-7)  for  p(v,x),  and  starting  values 


E(is,x)  =  2  #(-(2x)li)  ,  E  (  1 ,  x )  =  exp(-x)  . 


We  also  have  closed  form 


n-1 


E(n,x)  =  exp(-x)  )  x^/k!  for  nil. 

k=0 


(B-14) 


(B-15) 


(B-16) 


(B-17) 
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APPENDIX  C.  MOMENTS  OF  NORMALIZER  OUTPUT  FOR  SIGNAL  ABSENT 

For  ease  in  determining  the  value  of  threshold  v  to  use  in 

(49),  in  order  to  realize  specified  false  alarm  probability  P_, 

r 

it  is  helpful  to  know  the  moments  of  the  normalizer  output  y/yQ 
when  signal  is  absent.  Here,  we  will  determine  the  n-th  moment 
of  ratio  y/Yq: 

A*n  =  <Y/Y0)  =  Y  Y0  ,  (C-l) 

where  we  have  used  the  statistical  independence  of  y  and  yQ.  By 
reference  to  (24),  we  find 


L- 1 


’n  =  f  dx  x  n  - -  exP(  x/2-?  =  1-n)  1  for  L  >  n+1  .  (C-2) 

i  2b  ( L— 1 ) !  2n  (L-l) ! 


In  a  similar  manner,  for  signal  absent  (replacing  L  by  K), 


vn  .  f  d*  *n  1  exp(-x/2)  (K-H-n)l  2 
Y  “J„  2K  ( K-l )  !  '  <K-D  ! 


n 


( C-3 ) 


The  first  two  moments  of  ratio  y/Yq  are  then  given  by 


K 


(K  +  1)  K 


=  TTr-J  for  L  >  2  ,  p2  =  (L  i  !  )■(■£'  2)  for  L  >  3  .  (C-4) 


The  variance  of  the  normalizer  output  is  then  given  by 


Var(y/y  )  =  ^  =  (L  +  K~  1)  K —  for  L  >  3  ,  (C-5) 

°  (L  -  1)Z  (L  -  2) 


when  signal  is  absent. 
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APPENDIX  D.  PROGRAM  FOR  NORMALIZER  SIGNAL-TO-NOISE  RATIO 


The  program  in  this  appendix  computes  the  false  alarm 
probability  by  means  of  (49)  and  the  detection  probability  by 
means  of  (36).  Since  exceedance  distribution  function  value  (36) 
is  independent  of  p,  a  guess  at  p  is  required,  preferably  one 
that  leads  to  a  rapidly  convergent  series.  Here,  this  choice  has 
been  made  by  trial  and  error  and  by  observing  and  noting  values 
for  p  that  were  advantageous  for  "nearby"  parameter  values.  For 
the  purpose  of  ease  of  programming  and  for  interpretation  of  the 
recursions  in  the  listing  below,  we  have  utilized  the  following 
notation: 


A(p) 

=  P  ap 

for  p  >  1  , 

( D—  1 ) 

F(p) 

-  *  *p 

for  p  >  0  , 

(D-2) 

-  P 

emk 

Bl(m,k)  =  1 

_  P 

(D-3) 

2Qmk 

Qmk  ' 

22mk  ' 

Cp(m,k)  =  Bl(m,k)  +  p  Be(m,k)  ,  Bp(m,k)  =  Bl(m,k)P  1  ,  (D-4) 

Dp  =  K  (|  -  l)  [l  -  §)P  .  ( D—  5 ) 


The  program  below  allows  for  arbitrary  values  of  the 
deterministic  signal-to-noise  ratios  (D^/N^;  however,  it  was 
utilized  here  only  for  the  special  case  of  all  equal  to 

zero.  Also,  the  program  is  written  such  that  lD^m/N0l  are  held 
fixed  while  Ej/Nq  is  varied  to  meet  the  specified  detection 
probability  PD- 
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10 

!  NORMAL  1 2ER ,  NUWC-NPT  TR 

10275,  21  JANUARY  1993 

20 

Pf=.01  ! 

FALSE  ALARM  PROBABILITY 

30 

Pd*. 9  ! 

DETECTION  PROBABILITY 

40 

K*2  ! 

NUMBER  OF  SIGNAL  PULSES 

50 

L*10  ! 

NUMBER  OF  NOISE-ONLY  BINS 

60 

M*2  ! 

NUMBER  OF  FADING  COMPONENTS 

70 

DATA  l.,l.  ! 

POWER  RATIOS  r<m>  FOR 

1 

<  = 

m 

<*  M 

80 

DATA  0.  , 0. , 0. , 0.  ! 

Dkm/No  FOR  1  <=  k  <=  K, 

1 

<  = 

m 

<  *  M 

90 

DATA  1. ,2.  ! 

TIMES  tk  (SEC)  FOR 

1 

<  = 

k 

<  =  K. 

100 

DATA  0. , 0.  ! 

FREQUENCIES  fk  (HZ)  FOR 

1 

<  = 

k 

<  =  K 

110 

Elno0=20.  ! 

E 1  /No  STARTING  VALUE  (*0E5S) 

120 

Elnol=Elno0*.  1  ! 

E 1 ^No  INCREMENT 

130 

Ef *1 . E-15  ! 

TOLERANCE  ON  Pf 

140 

Ed*l . E-10  ! 

TOLERANCE  ON  Pd 

150 

P=500  ! 

MAXIMUM  NUMBER  OF  TERMS 

IN 

SUM 

(36) 

160 

OPTION  BASE  1 

170 

COM  DOUBLE  K,L,M,P  ! 

INTEGERS  (NOT  DOUBLE  PRECI 

SION) 

180 

DOUBLE  Ms, Ks, Js, LI  ! 

INTEGERS 

190 

DIM  Rs<5>,Dn<10,5),Ts<10),Fs<10),Psi (5>,Cbar(10, 10) 

200 

DIM  U<10,10),V<16,10),Eig<10),Sq<10,5> 

210 

COM  V,B2,Prod(5, 10),Es(5 

, 10) , E(0: 500) , B1 (5, 10) 

220 

COM  Be<5, 10), Bp <5, 10),Cp<5, 10) , F<0: 500) , AC 500) 

230 

REDIM  Rs<M),Dn(K,M>,Ts(IO,Fs(K),Psi (M) , Cbar(K, K) 

240 

REDIM  U(K,K),V(K,K),Eig<K),Sq(K,M) 

250 

REDIM  Prod<M,K),Es<M,K>, 

E(0:P),B1(M,K) 

260 

REDIM  Be<M,K),Bp(M,K),Cp(M,K),F<0:P>,A(P) 

270 

PRINT 

280 

PRINT  "Pf  * " ; P f ; "  K 

*  “  ;  K ;  "  L  =  "  }  L ;  "  M 

= 

H  • 

9 

M 

290 

READ  Rs<*),Dn<*),Ts<*),Fs<*)  !  Dn(*)  WAS  FILLED 

IN 

THE  ORDER: 

308 

S=0.  !  Dn< 1 , 1 ) 

, Dn( 1 , 2) , Dn( 1 , 3) . Dn(k 

,m), 

•  •  • 

, Dn(K, M) 

310 

FOR  Ms* 1  TO  M 

320 

S*S+Rs(Ms) 

330 

NEXT  Ms 

340 

FOR  Ms* 1  TO  M 

350 

Psi <Ms)=Rs(Ms)/S  ! 

Clj  (191)3 

360 

NEXT  Ms 

370 

FOR  Ks* 1  TO  K 

380 

FOR  Js*l  TO  K 

390 

Cov=FNCov(Ts(Ks)-Ts( Js) , 

Fs(Ks)-Fs< Js)  ) 

400 

Cbar (Ks, Js)*Cov  ! 

NORMALIZED  COVARIANCE  MATRIX 

1 1 

;  (127)3 

410  NEXT  J* 

420  NEXT  Ks 

430  MAT  U=Cbar 

440  CALL  Svd(  K ,K,U(*),V(*),Eig(*)) 

430  PRINT  "EIGENVALUES:" 

460  PRINT  Eig<*> 

470  FOR  Ms*  1  TO  M 

480  T*2. *Psi (Ms) 

490  FOR  Ks* 1  TO  K 

500  Prod<Ms , Ks > =  T*E i g< Ks ) 

510  Sq<Ks, Ms)*SQR(2. *Dn(Ks, Ms) ) 

520  NEXT  Ks 

530  NEXT  Ms 

540  FOR  Ks* 1  TO  K 

550  FOR  Ms* 1  TO  M 

560  S=0. 


!  OUTPUTS:  U<*),V(*),Eig<*) 
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570 
588 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850  fls 
860 
870 
830 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1  120 
1130 
1140 


FOR  Js=l  TO  K 

S=S+V< Js,Ks)*Sq< Js,Ms> 

NEXT  Js 
Es  <  Ms , Ks ) =S*S 
NEXT  Ms 
NEXT  Ks 

IF  L< 3  THEN  680 

Mean=K/(L-l )  !  (C-4) 

Var=(L+K-l  >*K/<  (L-l  W*(L-2>>  !  (C-5) 

V0=Mean+SQR<Var )*3.  !  THRESHOLD  STARTING  VALUE  < GUESS) 

GOTO  690 

V0=200  !  THRESHOLD  STARTING  VALUE  (ROESS) 

V 1 =  V  0  * .  1  !  THRESHOLD  INCREMENT 

CALL  Inversfunctionl<-Pf,Ef,V0,Vl,V> 

Beta=23.  !  SCALING  FACTOR  (GUESS) 

PRINT  "Threshold  V  =";V;"  Beta  =";Beta 

B2=Beta/2. 

X=V'(B2+V> 

Tk  =  (B2/'(B2  +  V))^L 
L1=L-1 

FOR  Ks= 1  TO  K-l 
Tk=Tk*X*(Ks+Ll)'Ks 
NEXT  Ks 

FOR  Ks=K  TO  K+P 
Tk=Tk*X*(Ks+Ll )/Ks 
E  <  Ks-K ) =Tk 
NEXT  Ks 
S  =  Tk 

FOR  Ks=K+P+ 1  TO  2000 

Tk=Tk*X*(Ks+Ll)/Ks 

S=S+Tk 

IF  Tk.<=S*l.E-15  THEN  920  !  ERROR  TOLERANCE 

NEXT  Ks 

PRINT  "2000  TERNS  ARE  INSUFFICIENT  IN  LINE  As" 

PAUSE 

E(P)=S 

FOR  Ks  =  P- 1  TO  0  STEP  -1 
E(Ks)=S=S+E(Ks> 

NEXT  Ks 

CALL  Inoersfunctior»2(Pd,Ed,Elno0,Elnol,Elro) 
Db=10.*LGT(Elno)  !  SNR  El^No  in  dB 

Y1=MIN(F<*>  > 

Y2=MAX(F<*> ) 

PRINT  "dB  =" ;Db;"  Y1  =  "jYlj"  Y2  =";Y2 

GINIT 

PLOTTER  IS  "GRAPHICS" 

GRAPHICS  ON 

WINDOW  0. ,50. , Y 1 , Y2 
LINE  TYPE  3 
GRID  10. , . 1 
LINE  TYPE  1 
FOR  Js=0  TO  P 

PLOT  J  s ,  F  <  J  s  )  !  EXPANSION  COEFFICIENTS  F  f(sub  p) 

NEXT  Js 

PENUP 

PAUSE 

END 

I 
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1150  DEF  FNCov  < T au , Nu )  !  NORMALIZED  COVARIANCE: 

1160  Rho= .  5  !  R1 1 <Tau, Nu)/Rl 1 <0, 0) ;  SEE  (4) 

1170  A=-.5*L0G<Rho) 

1180  B  =  0 . 

1190  Cov=EXP<-A*ABS<Tau)-B*ABS<Nu) >  !  EXPONENTIAL 

1200  RETURN  Coo 

1210  FNEND 

1220  ! 

1230  DEF  FNPf <V>  !  PROBABILITY  OF  FALSE  ALARM 

1240  COM  DOUBLE  K , L  !  INTEGERS 

1250  DOUBLE  Ks,Ll  !  INTEGERS 

1260  IF  V<=0.  THEN  RETURN  1. 

1270  V1=V/<1.+V> 

1280  Pf=T=EXP<-L*LOG<l.+V)) 

1290  L 1 =L- 1 

1300  FOR  Ks=  1  TO  K-l 

1310  T=T*V1*<L1+Ks)/Ks 

1320  Pf=Pf+T 

1330  NEXT  Ks 

1340  RETURN  -Pf  !  -  TO  YIELD  INCREASING  FUNCTION  OF  V 

1350  FNEND 

1360  ! 

1370  DEF  FNPd<Elno)  !  PROBABILITY  OF  DETECTION 

1380  COM  DOUBLE  K,L,M,P  !  INTEGERS 

1390  COM  V,B2,Prod<*),Es<*),E<*),Bl<*>,Be<*>,Bp(*>,CpC*>,F<*>,A(*) 

1400  DOUBLE  Ms , Ks , Ps  !  INTEGERS 

1410  To  1 = 1 . E-8  !  RELATIVE  ERROR  OF  SUM 

1420  F  =  1 . 

1430  S=0. 

1440  FOR  Ms*  1  TO  M 

1450  FOR  Ks*  1  TO  K 

1460  Q*1 . +Elno*Prod<Ms, Ks) 

1470  Eq=Es<Ms,Ks)/Q 

1480  F=F*Q 

1490  S=S+Eq 

1500  Bq=B2/Q 

1510  Bl<Ms,Ks)=l.-Bq 

1520  Be<Ms,Ks)=Bq*Eq 

1530  BpCMs, Ks)=l . 

1540  Cp  <  Ms , Ks )  =  1 . -Bq 

1550  NEXT  Ks 

1560  NEXT  Ms 

1570  F<0)=EXP(K*LOG<B2)-.5*S)/SQR<F) 

1580  Pd=l . -F<0)*E<0) 

1590  B2 1 = 1 . -B2 

1600  Dp=K*<M/2. -1 . ) 

1610  T 1  =  1 . 

1620  B:  FOR  PsM  TO  P 

1630  Dp=Dp*B2 1 

1640  S=0. 

1650  FOR  Ms* 1  TO  M 

1660  FOR  Ks* 1  TO  K 

1670  IF  P S *  1  THEN  1690 

1680  Bp<Ms,Ks)=Tl=Bp<'\s,Ks>*Bl<Ms,Ks> 

1690  Cp<Ms,Ks)=T2=Cp(Ms,Ks>+Be<Ms,Ks> 

1700  S=S+T1*T2 

1710  NEXT  Ks 

1720  NEXT  Ms 
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1730  A<Ps)=.5*S-Bp 

1740  8=0. 

1750  FOR  Ms=  1  TO  Ps 

1760  S=S  +  fl<Ms)*F<'Ps-Ms> 

1770  NEXT  Ms 

1780  F<Ps)=F  =  S/Ps  !  APPENDIX  FI 

1781  Do  =  De 1 

1790  De 1 =F*E  <  Ps ) 

1800  Pd=Pd-De 1 

1801  IF  Ps=l  THEN  1320 

1810  IF  ABS<Do)+ABS<Del  K=Pd*Tol  THEN  1850 

1820  NEXT  Ps 

1830  PRINT  P; "TERMS  ARE  INSUFFICIENT  IN  LINE  B:" 

1840  PAUSE 

1850  FOR  Ms=Ps+ 1  TO  P 

1360  F  <  Ms ) =A  <  Ms ) =0 . 

1870  NEXT  Ms 

1380  PRINT  "Pd  =" ; Pd; "  Ps  =";Ps, "El/No  =";Elno 

1890  RETURN  Pd 

1900  FNEND 

1910  i 

1920  SUB  I nversf  unc  t  i  on  1  ( Des  i  red ,  Err or ,  X 1 ,  De  1  ,  X2  > 

1930  X2=X 1  +De 1 

1940  F 1 =FNPf <  X 1  )  !  FALSE  ALARM  PROBABILITY 

1950  F2  =  FNPf  <  X2 ) 

1960  IF  F2  >  =  Des i red  THEN  2010 

1970  X 1 =  X  2 

1930  X2=X2+De 1 

1990  F 1 =F2 

2000  GOTO  1950 

2010  IF  F 1 <  =Des i red  THEN  2070 

2028  X2=X 1 

2030  Xl=Xl-Del 

2040  F2=F 1 

2050  F l=FNPf <X1 ) 

2060  GOTO  2010 

2070  Xa=Xl 

2080  Xb=X2 

2090  IF  (.  F2-Desi  red)<  CDs  si  red -FI)  THEN  2160 

2100  T=X 1 

2110  X 1 =X2 

2120  X2=T 

2130  T=F 1 

2140  F 1 =F2 

2150  F2=T 

2160  IF  ABS < F2-Des i r ed > < Err  or  THEN  2260 

2170  IF  F2  =  F  1  THEN  2260 

2130  T=<X1*<F2-Desi red)-X2*<Fl -Desired) )/<F2-Fl ) 

2190  T  =  MAX  <  T , Xa) 

2200  T=MIN<T, Xb> 

2210  X 1 =X2 

2220  X2=T 

2230  F 1 =F2 

2240  F2=FNPf  <  X2 ) 

2250  GOTO  2160 

2260  SUBEND 

2270  i 
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2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
28  30 
2340 
2350 
2860 


SUB  Inversfunction2<Desired,Error,Xl,Del,X2) 
X2=X 1 +De  1 

Fl=FNPd<Xl )  !  DETECTION  PROBABILITY 

F2=FNPd<X2> 

IF  F2  >=Des i red  THEN  2370 

X 1  =X2 

X2=X2+De 1 

F1=F2 

GOTO  2310 

IF  F l<=Desi red  THEN  2430 
X2  =  X  1 
X 1  =X  1 -De 1 
F2  =  F  1 

F l=FNPd<Xl ) 

GOTO  2370 

Xa=Xl 

Xb=X2 

IF  <F2-DesiredX<DesireJ-Fl>  THEN  2520 
T  =  X1 
X 1  =X2 
X2  =  T 
T  =  F  1 
F 1  =F2 
F2  =  T 

IF  ABS<F2-DesiredXError  THEN  2620 
IF  F2=F 1  THEN  2620 

T=<X1*<F2-Desi red)-X2*<Fl-Desi red) )/<F2-Fl ) 

T  =  MAX  <  T , Xa) 

T=MIN<T,Xb) 

X1=X2 
X2  =  T 
F 1  =F2 

F2=FNPd<X2) 

GOTO  2520 
SUBEND 


SUB  Svd<  DOUBLE  M , N, REAL  A  <  *  )  ,  V  <  *  )  ,  W  <  *  )  ) 

THIS  SUBROUTINE  COMPUTES  THE  SINGULAR  VALUE  DECOMPOSITION 
OF  AN  ARBITRARY  REAL  MxN  MATRIX  A;  A  =  U  W  Vt,  M  >=  N. 


U  IS  MxN,  V  IS  NxN ,  U  IS  NxN : 
ALLOCATE  Rvl<l:N> 

IF  MXN  THEN  2720 

PRINT  " M < N  IS  DISALLOWED" 

PAUSE 

DOUBLE  I, J,K,L, Its,Nm, Jj 
G=Sc  al e=Anor m=0 . 

FOR  1=1  TO  N 
L=  I  + 1 

Rwl < I )=Scal e*G 
G=S=Sc al e=0 . 

IF  I >M  THEN  3060 
FOR  K= I  TO  M 
Scale=Scal*+ABS<A<K, I>) 

NEXT  K 

IF  Sc  al e  =  0 .  THEN  3060 

FOR  K= I  TO  M 

Aa=A^K,,  I  )  =  A < K ,  I  >  -  Seal  e 

S=S+Aa* Aa 

NEXT  K 


=  DIAG<D<n)>. 

NUMERICAL  RECIPES,  PAGES  60-64 
A<*>  IS  OVER-WRITTEN 


INTEGERS  < NOT  DOUBLE  PRECISION) 
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2870  F=fi<I,I> 

2880  G=-SQR  <  S  > 

2890  IF  F<0.  THEN  G=-G 

2900  H=F*G-S 

2910  fl< I , I )=F-G 

2920  IF  I =  N  THEN  3030 

2930  FOR  J=L  TO  N 

2940  S=0. 

2950  FOR  K= I  TO  M 

2960  S=S+fl<K, I >*A<K, J) 

2970  NEXT  K 

2980  F=S^H 

299b  FOR  K= I  TO  M 

3000  fi<K, J>+F*fi<K, I> 

3010  NEXT  K 

3020  NEXT  J 

3030  FOR  K=  I  TO  (1 

3040  fKK,  I>=A<K,  I)*Scale 

3050  NEXT  K 

3060  W< I )=Scal e*G 

3070  G=S=Sc  al e=0. 

3080  IF  COM)  OR  <  I  =N )  THEN  3380 

3090  FOR  K=L  TO  N 

3100  Seal e=Sca1 e  +  fiBSCfi< I ,  K> ) 

3110  NEXT  K 

3120  IF  Seal e=0.  THEN  3380 

3130  FOR  K=L  TO  N 

3140  fia=fl<I,K)=fl<I,K)/'Sca1e 

3150  S=S+fla»fla 

3160  NEXT  K 

3170  F=fl< I , L) 

3180  G=-SQR<  S ) 

3190  IF  F<0.  THEN  G=-G 

3200  H=F*G-S 

3210  fi< I , L)=F-G 

3220  FOR  K=L  TO  N 

3230  R«v>l<K)=fl<I,K)/'H 

3240  NEXT  K 

3250  IF  I *M  THEN  3350 

3260  FOR  J=L  TO  M 

3270  S=0. 

3280  FOR  K=L  TO  N 

3290  S  =  S  +  PKJ,K>*fl<I,K) 

3300  NEXT  K 

3310  FOR  K=L  TO  N 

3320  fl<J,K)=fi<J,K)+S*Rvl<K) 

3330  NEXT  K 

3340  NEXT  J 

3350  FOR  K=L  TO  N 

3360  fl<I,K>=A<I,K>*Scale 

3370  NEXT  K 

3380  flnorm*MAX<Anorm, ABS<W< I ) )+ABS<Rvl < I ) ) ) 

3390  NEXT  I 

3400  FOR  I =N  TO  1  STEP  -1 

3410  IF  I >=N  THEN  3580 

3420  IF  G*0 .  THEN  3550 
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3430  FOR  J=L  TO  N 

3440  V<J, I)=fl(I, 

3450  NEXT  J 

3460  FOR  J  =  L  TO  N 

3470  S  =  0. 

3480  FOR  K=L  TO  N 

3490  SaS+fl<I,K)*V<K,J) 

3500  NEXT  K 

3510  FOR  K»L  TO  N 

3520  VCK, J)=V<K, J>+S*V<K, I> 

3530  NEXT  K 

3540  NEXT  J 

3550  FOR  J=L  TO  N 

3560  V<I , J)=V< J, I )=6. 

3570  NEXT  J 

3580  V  < I , I >  =  1 . 

3590  G=Rv 1 < I ) 

3600  L-I 

3610  NEXT  I 

3620  FOR  I =N  TO  1  STEP  -1 
3630  L* I + 1 

3640  G*W  Cl) 

3650  IF  I >=N  THEN  3690 

3660  FOR  J*L  TO  N 

3670  fl<I,J)=0. 

3680  NEXT  J 

3690  IF  G=0 .  THEN  3860 

3700  G=1./'G 

3710  IF  I =N  THEN  3820 

3720  FOR  J=l  TO  N 

3730  S*0. 

3740  FOR  K=L  TO  M 

3750  S=S+fi<K, I J) 

3760  NEXT  K 

3770  F=S/R< I , I >*G 

3780  FOR  K= I  TO  M 

3790  R<K,J)=fl<K,J)+F*R<K, I> 

3800  NEXT  K 

3810  NEXT  J 

3820  FOR  J=I  TO  M 

3830  fl< J, I)*fl< J, I )*G 

3840  NEXT  J 

3850  GOTO  3890 

3860  FOR  J=I  TO  tl 

3870  fl<J,I)=0. 

3880  NEXT  J 

3890  fl<I,I)*fl<I,I>+l. 

3900  NEXT  I 

3910  FOR  K=N  TO  1  STEP  -1 

3920  FOR  1 1  s =  1  TO  30 

3930  FOR  L=K  TO  1  STEP  -1 

3940  Nm=L- 1 

3950  IF  <RBS<Rvl <L) )+flnorm)=flnorm  THEN  4170 

3960  IF  <RBS(W<Nm) )+Pnorm)=Rnorm  THEN  3980 

3970  NEXT  L 

3930  C*0. 

3990  S* 1 . 
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4000  FOR  I=L  TO  K 

4010  F=S*Rv 1 < I ) 

4020  Rvl < I )  =  C*Rul<  I  ) 

4030  IF  < ABS < F ) +flnorm ) =Anorm  THEN  4170 

4040  G=U< I > 

4050  H=SQR<F*F+G*G) 

4060  W< I )=H 

4070  H=l./-H 

4080  C=G*H 

4090  S=-F*H 

4100  FOR  J=1  TO  II 

4110  Y-A< J, Nm) 

4120  Z*R ( J ,  I ) 

4130  A<J,Nm)=Y*C+Z*S 

4140  A< J, I)=-Y*S+Z*C 

4150  NEXT  J 

4160  NEXT  I 

4170  Z=W<IO 

4180  IF  LOK  THEN  4250 

4190  IF  Z>=0.  THEN  4240 

4200  W<K)=-Z 

4210  FOR  J=1  TO  N 

4220  V<J,K)=-V<  J,K). 

4230  NEXT  J 

4240  GOTO  4780 

4250  IF  1 1 s< 30  THEN  4280 

4260  PRINT  "NO  CONVERGENCE  IM  30  ITERATIONS" 
4270  PAUSE 

4280  X  =  W  <  L ) 

4290  Nm=K-l 

4300  Y=W<Nm) 

4310  G=Rvl <Nm) 

4320  H  =  Rv 1  <  K) 

4330  F=<(Y-Z)*<Y+Z)+<G-H)*<G+H) )/<2. *H*Y) 

4340  G=SQR<F*F+1. > 

4350  Aa=ABS<G> 

4360  IF  F<0.  THEN  Aa=-Aa 

4370  F=<<X-Z)*<X  +  Z)+H*<  <Y/<F  +  Aa))-H>)^X 

4380  C=S= 1 . 

4390  FOR  J  =  L  TO  Nm 

4400  I=J+1 

4410  G*Rv 1 < I ) 

4420  Y=W<I> 

4430  H=S*G 

4440  G=C*G 

4450  Z=SQR<F*F+H*H) 

4460  Rv 1 <  J ) =Z 

4470  C*F/”Z 

4480  S=H/Z 

4490  F=X*C+G*S 

4500  G=-X*S+G*C 

4510  H=Y*S 

4520  Y=Y*C 

4530  FOR  Jj*l  TO  N 

4540  X*V<Jj,J> 

4550  Z=V<Jj,I) 

4560  V< Jj , J)»X*C+Z*S 

4570  V< Jj, I )=-X*S+Z*C 

4580  NEXT  Jj 
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4590  Z=SQR<F*F+H*H) 

4600  W< J)=2 

4610  IF  2=0.  THEN  4650 

4620  2=1. /Z 

4630  C=F*2 

4640  S  =  H*Z 

4650  F  =C*G+S*Y 

4660  X=-S*G+C*Y 

4670  FOR  Jj=l  TO  M 

4680  Y=fl(Jj,J) 

4690  2=fl(Jj,I> 

4700  fl< Jj , J>=Y*C+Z*S 

4710  fl< Jj , I )=-Y*S+Z*C 

4720  NEXT  Jj 

4730  NEXT  J 

4740  Rv 1 <  L ) =0 . 

4750  Rvl <K>=F 

4760  W<  K  >  =X 

4770  NEXT  Its 

4780  NEXT  K 

4790  SUBEND 


The 

values  of  threshold 

v  and  required 

input 

signal -to-noise 

ratio  E 

./N  ,  that  are 
1  o 

:  yielded  by  the 

above 

program, 

are  listed 

below. 

as  a  function 

of 

L, 

the  number 

of  noise-only 

bins.  The 

values 

of  scaling  factor 

fi 

are  chosen 

to  approximately  minimize 

the 

number  of  terms. 

V 

used  in  the 

p-summation 

in 

(36)  . 

L 

V 

E1/Nq  (dB) 

£ 

NP 

1 

198.5 

27.263 

535 

14 

2 

15.98 

19.123 

85 

12 

3 

6.099 

16.566 

48 

10 

4 

3.503 

15.317 

36 

10 

5 

2.398 

14.577 

31 

10 

6 

1.804 

14.086 

28 

10 

7 

1.439 

13.738 

26 

10 

8 

1.193 

13.477 

25 

9 

9 

1.018 

13.274 

24 

9 

10 

.8861 

13.113 

23 

9 

16 

.4962 

12.569 

21 

8 

24 

.3115 

12.267 

20 

8 

32 

.2267 

12.117 

19 

8 

48 

.  1467 

11.966 

18 

8 

64 

.  1084 

11.891 

18 

8 

100 

.0683 

11.810 

18 

8 
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